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Abstract 

According to extensive experimental findings, the Ginzburg temperature to for ionic fluids differs 
substantially from that of nonionic fluids [Schroer W., Weigartner H. 2004 Pure Appl. Chem. 76 

O ' 

19] . A theoretical investigation of this outcome is proposed here by a mean field analysis of the 

^— > ■ 

interplay of short and long range interactions on the value of tQ. We consider a quite general 
continuous charge-asymmetric model made of charged hard spheres with additional short-range 
interactions (without electrostatic interactions the model belongs to the same universality class 

<-> ■ 

as the 3D Ising model). The effective Landau-Ginzburg Hamiltonian of the full system near its 
gas-liquid critical point is derived from which the Ginzburg temperature is calculated as a function 
of the ionicity. The results obtained in this way for tc are in good qualitative and sufficient 

in 

O ■ quantitative agreement with available experimental data 
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I. INTRODUCTION 



It is known that electrostatic forces determine the properties of various systems: phys- 
ical as well as chemical or biological. In particular, the Coulomb interactions are of great 
importance when dealing with ionic fluids i.e., fluids consisting of dissociated cations and 
ions. In most cases the Coulomb interaction is the dominant interaction and due to its 
long-range character can substantially affect the critical properties and the phase behavior 
of ionic systems. Thus, the investigations concerning these issues are of great fundamental 
interest and practical importance. 

Over the last ten years, both the phase diagrams and the critical behavior of ionic solu- 
tions have been intensively studied using both experimental and theoretical methods. These 
studies were stimulated by controversial experimental results, demonstrating the three types 
of the critical behavior in electrolytes solutions: (i) classical (or mean-field) and (ii) Ising-like 

nnnrtnn 

behavior as well as (iii) crossover between the two [1|, [2|, [3|, |4|, [5|, [6j . In accordance with these 
peculiarities, ionic solutions were conventionally divided into two classes, namely: "solvo- 
phobic" systems with Ising-like critical behavior in which Coulomb forces are not supposed 
to play a major role (the solvent is generally characterized by high dielectric constant) and 
"Coulombic" systems in which the phase separation is primarily driven by Coulomb in- 
teractions (the solvent is characterized by low dielectric constant). Hence the criticality 
of the Coulombic systems became a challenge for theory and experiment. A theoretical 
model which demonstrates the phase separation driven exclusively by Coulombic forces is 
a restricted primitive model (RPM) [7J, |8j. In this model the ionic fluid is described as 
an electroneutral binary mixture of charged hard spheres of equal diameter immersed in a 
structureless dielectric continuum. Early studies ^, Q, 11 1 established that the model has a 



gas-liquid phase transition. A reasonable theoretical description of the critical point in the 



RPM was accomplished at a mean-field (MF) level using integral equation methods |8j, [12 
and Debye-Huckel theory 13j. Due to controversial experimental findings, the critical be- 
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havior of the RPM has been under active debates 
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strong evidence 
and theoretical 



'or an Isin g un iversal class has been found by recent simulations 
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29j studies. 



In spite of significant progress in this field, the criticality of ionic systems are far from 
being completely understood. The investigation of more complex models is very important 
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in understanding the nature of critical behavior of real ionic fluids demonstrating both the 
charge and size asymmetry as well as other complexities such as short-range attraction. A 
description of a crossover region when the critical point is approached is of particular interest 
for such models. Based on the experimental findings one can suggest that in ionic fluids 
the temperature interval of crossover regime, characterized by the Ginzburg temperature, 
is much smaller than observed in nonionic systems 



reported for the systems Na — NH3 



301 ] (see also 



31 



5|. In particular, a sharp crossover was 



23, 



33(). The analysis of experimental 



data for various ionic solutions confirmed that such systems generally exhibit crossover 
or, at least a tendency to crossover from the Ising behavior asymptotically close to the 
critical point, to the mean-field behavior upon increasing distance from the critical point 
34J] . Moreover, the systematic experimental investigations of the ionic systems such as tetra- 
n-butylammonium picrate, Bu4NPic, (for tetra-n-butylammonium picrate we will follow the 
notations from 5|, |6j) in long chain n-alkanols with dielectric constant ranging from 3.6 for 
1-tetradecanol to 16.8 for 2-propanol suggest an increasing tendency for crossover to the 
mean-field behavior when the Coulomb contribution becomes essential 3, H, [35I . They also 
indicate that the "Coulomb limit" reduced temperature of the RPM TL— 0.05 is valid for 



the almost non-polar long chain alkanols 



351 ] . It has been stressed |35| that for solutions 



of Bu4NPic in 1-alkanols, the upper critical solution points are found to increase linearly 
with the chain length of the alcohols (that corresponds to the decrease of dielectric constant 
of the solvent). The experimental data for the critical points and the dielectric permittivities 
for solutions of Bu4NPic in 1-alkanols are given in Table 1 35 1. 

Theoretically the crossover behavior in ionic systems was firstly studied for the RPM 
HQ . Il6j]. The results obtained for the Ginzburg temperature were similar to those found 
for simple fluids in comparable fashion that is in variance to what is expected from the ex- 
periments 5|, |6j. Nearly at the same time in 36j the crossover behavior of the lattice version 
of a fluid exhibiting the Ising behavior was studied as additional symmetrical electrostatic 
interactions were turned on. Based on the microscopic ground, the effective Hamiltonian 
in terms of the fluctuating field conjugate to the number density was derived in this work. 
Then, the crossover between the mean-field and Ising-like behavior was estimated using the 
Ginzburg criterion. The resulting crossover temperature calculated as function of the ionicity 
X, which defines the strength of the Coulomb interaction relative to the short-range inter- 
action, indicates its weak dependence but with the trends correlating with those observed 
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TABLE I: The experimental parameters of the critical points (critical temperature T c , critical mass 
fraction w c ) and the corresponding dielectric constants e for solutions of Bu4NPic in 1-alkanols 351 ] . 



Solvent 


<T C ) 


TJK 


Wc 


1-oktanol 


9.5 


298.55 


0.336 


1-nonanol 


7.9 


308.64 


0.325 


1-decanol 


6.4 


318.29 


0.3152 


1-undecanol 


5.4 


326.98 


0.303 


1-dodecanol 


4.7 


335.91 


0.2951 


1-tridecanol 


4.3 


342.35 


0.284 


1-tetradecanol 


3.6 


351.09 


0.2721 



experimentally. 

In this paper we are also interested in the critical behavior of ionic fluids. In particular, we 
study the effect of the interplay of short-range and long-range interactions on the crossover 
behavior in such systems. We consider a continuous version of the charge-asymmetric ionic 
fluid in which both the long-range Coulomb and short-range van-der Waals-like interactions 
are included. Following [36J we introduce the ionicity 

J= 1 = ^L (1) 

where qi is the charge on ion i, kg is the Boltzmann constant, T is the temperature, a is 
collision diameter and e is the dielectric constant. Then we derive the effective Hamiltonian 
of the charge-asymmetric model in the vicinity of the gas-liquid critical point. As in 36], 
the coefficients obtained for the effective Hamiltoninan have the forms of expansions in the 
ionicity but with new terms that appear in this case. Based on this Hamiltonian we estimate 
the Ginzburg temperatures as functions of the ionicity. 

The layout of the paper is as follows. In Section 2 we introduce a continuous charge- 
asymmetric model with additional short-range attractive interactions included. We derive 
here the functional representation of the grand partition function of the model in terms of the 
fluctuating fields </?f and (p£ conjugate to the total density and charge density, respectively. 
Section 3 is devoted to the derivation of the effective GLW Hamiltonian in the vicinity of 
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the critical point. In Section 4 we calculate the Ginzburg temperature as a function of the 
ionicity for different values of the range of the attractive potential. We conclude in Section 5. 



II. BACKGROUND 
A. Model 

Let us start with a general case of a classical charge-asymmetric two-component system 
consisting of N particles among which there are Ni particles of species 1 and N2 particles 
of species 2. The pair interaction potential is assumed to be of the following form: 

M^CW + ffl + tt). (2) 

where (r) is the interaction potential between the two additive hard spheres of diameters 
a a and a p. We call the two-component hard sphere system a reference system. Thermody- 
namic and structural properties of the reference system are assumed to be known. 0ff (r) 
is the potential of the short-range (van-der-Waals-like ) attraction. 0?g(r) is the Coulomb 
potential: (p^pij-) = q a q/3(p C (r)/e, where c (r) = 1/r and e is the dielectric constant. The 
solution is made of both positive and negative ions so that the electroneutrality condition 
is satisfied, i.e. Y^a=i Qa c ® = 0> where c a is the concentration of the species a, c a = N a /N. 
The ions of the species a — 1 are characterized by their hard sphere diameter a 1 and their 
electrostatic charge +q and those of species a — 2, characterized by diameter <t 2 , bear 
opposite charge — zq§ (qo is elementary charge and z is the parameter of charge asymmetry). 
In general, the two-component system of hard spheres interacting via the potential 
can exhibit both the gas-liquid and demixion critical points which belong to the 3D Ising 
model universal class. 

We consider the grand partition function (GPF) of the system which can be written as 
follows: 



= L L 11 AM / (dr)eXp 

A r i>0 7V 2 >0a=l,2 a ' 



a/3 ij 



(3) 



Here the following notations are used: u a is the dimensionless chemical potential, u a = 
fi\i a — 3 In A, fi a is the chemical potential of the ctth species, (5 is the reciprocal temperature, 
A -1 = (2nm a f3~ l /ti 2 ) 1 / 2 is the inverse de Broglie thermal wavelength; (dr) is the element 
of configurational space of the particles: (dr) = Yl a dl", dr a = dr^dr^ . . . dr£j a . 



Let us introduce the operators and p^ 



pI 



QaPk.,a: 



which are combinations of the Fourier transforms of the microscopic number density of the 
species a: p^ a = ^ exp(— ikr"). In this case the part of the Boltzmann factor entering eq. 
([3]) which does not include hard sphere interactions can be presented as follows: 



exp 







a/3 i,j 



exp 



2V 



+tnB(k)p°p D k + 2$ 3D {k)pip° h ) +4-J2 N « £(&) + 



where 



4>ss(k) 

4>DD{k) 

4>SD(k) 



1 

1 



z 2 ^ R {k) + 2z$ R (k) + ~^{k) 

4>li{k) ~ 2$?{k) + 4> S 22 R (k)] + 4>°(k) 
z4> s n R (k) + (1 - z)tf R (k) - $*{kj 



(4) 



(5) 



(1 + ^) 2 

with 4>*a'(k) being a Fourier transform of the corresponding interaction potential defined by 



fcS'(b)= / dr0^-(r)exp(-ikr) 



Now we simplify our model assuming that 

• The hard spheres will all be of the same diameter a a = a. 

• = = j%*(k) = SR (k). 

With these restrictions the uncharged system can only exhibit a gas-liquid critical point and 
a possible demixion is ruled out. 

Taking into account the assumptions mentioned above we thus have 



(f>ss(k) = 4> bK {k) < 0, <j> DD {k) = 4>°{k) > 0, (j)sD(k) = 0. 

Finally it will be convenient to introduce the effective range bsR of short-range interactions 
through the relations 



lSR 



(k) 



lSR 



(0) (1 - (b SR kf) + O(k') . 



(6) 
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B. Functional representation of the grand partition function of an ionic model 

Let us take advantage of the properties of Gaussian functional integrals to rewrite 



with 



5 

k 



N ws = J (dip 8 ) exp \-± ^ [ w s(k)\ 1 



and 



(d^) = n d ^ = n d^d^), A— s,d. 

k k 

In the above equations we also introduced the notations ws(k) = —P<j>ss(k)/V and 
w c (k) = 4> c {k)/V. 

As a result, we can rewrite S[i/ Q ] in the form of a functional integral 

= V7 tt" / (d^)(dOexp(-W[^,^ 5 ,^]) , (7) 
where the action Ti reads as 



k k 

-ln-HslVs + ^^ + ip 1 / 2 ^}, (8) 



where the " renormalized" chemical potentials v a are defined as 

^ = "« + ^E (-™s(k) +Pq 2 a w c (k)) , a = 1,2. (10) 

k 

Let us define Au s = V s — ^q and = Az/ S ' + (p^ with ip^ chosen as the chemical potential 
of the hard spheres. This leads to the relation 

u s + (p S = Vo+fi S - (11) 
Now we present lnH#s[. . .] in the form of a cumulant expansion 

lnE HS [...} = E^E E ^Vo^h,...,^...^ 

n>0 i„>l ki,...,k n 

•••^k n 4 1+ ...+k„, (12) 



l n + l 

where 97t^ [yoj^-Dj • • • , k n ] is the nth cumulant (or the nth order truncated correlation 
function) defined by 

Qn 1 _ c I" j 

TO!-. W, FDi = g- £ g-gjgg- ^ (13) 

In particular it follows from (TT5]) that 

OJtf =lBS ff5 [^,^]- (14) 

The expressions for the cumulants of higher order (for i n < 4) are given in Appendix A. 
It should be noted that, contrary to [36], ( TT2l) includes all powers (even and odd) of the 
field ip^ conjugate to the total number density. It should be clear that the coefficients in the 
cumulant expansion (12) depend on the chemical potential (or, equivalently, on the density). 

III. EFFECTIVE HAMILTONIAN IN THE VICINITY OF THE CRITICAL 
POINT 

Taking into account ffl2|) we can rewrite ©-([8]) as follows 

= -S-exp(-7?) ^(d^expU^^rVk^-k 

n>l ki,...,k n 

x5 kl +... +k jV[^], (15) 



S 



where 

n= 1 -[w s (o)r 1 (^ s ) 2 -^Hs[^}, 



V[^] = ^J{^ D )e W (-\^[w c {kT^ 

ki,k2,k3 ki,...,k4 

x$LS^+...+u + l E ^i 3) ^k 1 ^k 2 ^k3^ 4 ^ 1+ ...+k 4 + ...) • (16) 

ki,...,k 4 / 



It is worth noting here that unlike to the case considered in 36j we obtain in (fTEj) terms 
proportional to (v? 5 ) 2 (</5 D ) 2 and tp s (<£> D ) 3 . While the former is connected with an absence 
of a lattice symmetry, the the latter stems from charge asymmetry. 

Our aim now is to derive the effective Landau- Ginzburg (LG) Hamiltonian. Since we are 
interested in the gas-liquid critical point, this Hamiltonian should be written in terms of 
fields ip^. conjugate to the fluctuation modes of the total number density. 

To this end we integrate out cp^ in (|TB|) using a Gaussian measure. As a result, we can 
present V[£>j[] as follows: 



V[5lp 



■S] _ 

J V W( -, 



(17) 



where (■ ■ -)g means 



(•••)g = ^ /(<!/') ...ex,, | --} :Wr(M^^ k 
with VFc(fc) given by 

W c (k) = [w c (k)}~ 1 + y 2 Gi (18) 

and ?/ 2 being the ionicity introduced by (pQ): y 2 =1 

Taking into account (CQ) and the recurrence formulas of Appendix A A may be written 
as a formal expansion in terms of y 

^ = -y E ^(fo, ^ + fc3)<^k 1 ¥ 3 k 2 ^k3 5 k 1 +k 2 +k 3 - "J G 3^1, &2, &3 + h) 

ki,k2,k 3 ki,...,k4 

X<^k 2 ^ 3 ^k 4 ^k 1 + ...+k 4 " ^f^T^ E ^^i' fca + k * + ^ 



,k 4 

J> J>,„D S 

<3- 



x««^k 1+ ...+k 4 + .... (19) 

9 



In (I18 p - (jl9p the "tilde" over was omitted for the sake of simplicity. 

It should be mentioned that the dependence of G n (ki, k 2 , ■ ■ ■ , k n ) on the ki is very compli- 
cated. Since we consider here the behavior of the system near the critical point the limiting 
case of ki = is of particular interest. Therefore, we substitute in (JTTJ) 



G n {fa, k 2 , . . ■ , k n ) = G n (0, . . .) n>3 



and 



2l.2\ 



with 



G 2 (k)=G 2 (0)(l + g 2 k 



n2 _G 22 (0) ~ ( ^_ d 2 G 2 (k) 



2G 2 (0)' dk 
Having integrated out eq. (fTTj) 5[i/ a ] takes the form 



-LlKl + yi^HsWcik))- 1 [ (dy> 5 ) exp (-ft e/ Vl) 



(20) 



(21) 



n 

where we have for the coefficients a r 



n! 

n>0 ki,...,k n 



(22) 



Q-0 = ~~ 



2 4 / 

ai = (iV) H5 + ^ 5 (0)]- 1 A^-|-G 2 (0)^A(g) + ^f3G3(0) + 



(23) 



'l-z) 2 -2z 



xG 2 (0) 



£ A (<?) 



(24) 



a 2 = -[^r^G^^-l-^W^A^ + ^^Wl^A^Adk + ql), (25) 



2 



a 3 = G 3 (0)-^G 4 (0)^A(g) + -y 4 G 2 (0)G 3 (0)^A 2 (g), 



(26) 



a 4 = G 4 (0)-^G 5 (0)^A(g) + ^ 4 (3[G 3 (0)] 2 + 4G 2 (0)G 4 (0))^A 2 (g), (27) 



and the propagator A(q) is written as 

A(q)=A(q;y 2 ) = [Wc(q)]- 1 



w c {q) 



l + y 2 (N) HS w c (q) 



(28) 
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Coefficients (I23l)- fl27j) have the form of a formal expansion in terms of the ionicity X = y 2 . 
In our study all terms which do not exceed the fourth order of y are kept. The ionicity is 
small enough for large values of the dielectric constant and increases with its decrease. From 
this point of view we can consider the expansions in (T2l?l) - f[2"Tl) for large values of y 2 only as 
formal ones. It should be also noted that A(g) ~ 1/y 2 (see ( 1281) ) for large values of y 2 . 
Let us introduce 

r S R = [G 2 (0)^(0)]- 1 - 1 = (29) 



c.O 



where T c = T C (I = 0) is the mean- field critical temperature of the uncharged system. 
Taking into account (1291) we can rewrite —7i e ^ as follows: 



n^[lf s ] = -^[^(Or^A^^lnS^^-l^^o + ro^ 2 )^ 



s 

-k 



^0 \ ^ _s —S —S r u \ ^ 



"o 



ki,k 2 ,k 3 



ki,k2,k 3 ,k4 



x5ki+k 2 +k 3 +k 4 — h <p$, 



(30) 



where = V ^2(0)^ an d the following notations were introduced: 



ho 



rsR + 



r 



Tn 



''0 



Uq = — 



2 G 2 (0) 



£A( g )-^G 2 (0)£> 2 ( g ) 



4-7G 2 (0)^A( 9 )A< 2 )( g ) 



G 3 (0) y 2 G 4 (0) 



^A(g) + ^G 3 (0)^A 2 (g) 



[G 2 ]V2 ^G 2 (0) 2 g 2 ( ) ^ 

G 2 \G 2 (0) 2 G 2 (0)^ 2 q ^(0) 



x 



3[G 3 (0)] 2 + 4G 2 (0)G 4 (0) 



1/2 



~[G 2 

(l-z) 2 -2z 



(N) HS + [w s (0)}- l A^ _ y^A(q) + V - 



G 2 (0) 



, G0 

G 2 (0) 



-1 2^ 



(31) 
(32) 

(33) 



(34) 



(35) 



where 



— T, 



0' + 



°SR 



' ' G 2 (0)[w s (0)} 2 



(36) 
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with b 2 SR = b 2 SR w s {0) and A( 2 )(g) = <9 2 A(| k + q |)/<% 2 | fc=0 . 
Finally, we present (130!) as follows: 

-U^m = -^[^(OF^A^^lnH^^-i^^ + rV)^ 

k 

V^/HS ki.k2.k3 ^ /H5 ^,...^4 

x5 kl+k2 +k 3 +k4 - H^hs^q, (37) 

with 



r = r , r =r , u = — (iV)^ u = — (N) HS , h = h (N) H < 



At the critical point the following equalities hold 

r = 0, v = 0, h = 0, 

which give the equations for the critical parameters i.e., the temperature, the density and 
the chemical potential at the critical point. 

Equation (1371) gives the effective GLW Hamiltonian of the system (|2j) in the vicinity of 
the critical point. We are now in position to extract from eq. ( 1371) the Ginzburg temperature 
as a function of the ionicity. 

Now let us specify the short-range attraction, (p SR (r), in the form of the square-well 
potential 

0, < r < a 

—£, o < r < Act . 

0, r > Act 

It is worth noting here that the system of hard spheres interacting through the potential 
<f) SR (r) with A = 1.4 — 1.7 reasonably models most simple fluids [3?]]. The Fourier transform 
of <fi (r) for the case of the Weeks- Chandler- Andersen (WCA) regularization inside the 



</> SR (r) 



hard core 



has the form: 



<p SR (k) = 5fi (O)-^— [-Ax cos(Ax) + sin(Ax)], (38) 
(Xx) 6 

where x = ko and 4> SR (0) = — £cr 3 ^ [ A 3 . 

To be consistent we also use the WCA regularization scheme for the Coulomb potential 
which yields 

(f) c (x) = 4irsm(x)/x 3 . (39) 
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IV. GINZBURG TEMPERATURE 



Following 36( we can present the Ginzburg temperature by 



^''''^ [ii^n) ' (40) 

but in our case all quantities w, to and r 2 should be estimated at critical density r) c (y): 

u(y) = u(r] c (y),y), t (y) = t (rj c (y) , y) , r 2 (y) = r 2 {iq c {y),y). 

The density 77 enters the expressions for u, t and r 2 through the structure factors S n . A 
well-known criterium by Ginzburg predicts that the mean-field theory is valid only when 
tc «\ t | where t = T ~^^ and T c (y) are the mean-field reduced temperature and the 
mean-field critical temperature of the charged system at r\ — rj c (y), respectively. 

In PUI) t (y) measures the increase of the mean-field temperature of the charged system 
in respect to the uncharged system 

to(y) = ^ - 1 (41) 

which, for the model under consideration has the form: 

^ = -y^TnT + P^)J2^- ( 42 ) 

Z G 2(U) q Z q 

Taking into account that 

i t 



G 2 (0)wc(0) 

and equation (I4T1) . we can rewrite ( J36l) as follows: 



^c,0 



6 : 



Tc 



/ + -%(l+to(y)). (43) 



5 * a ' wc(0) 
For the uncharged model the Ginzburg temperature reduces to 

1 ^ 2 
ta(X=0))- 4 



^ Sl[r SR (l = 0)f 

where S n is given by (USD and r SR (l = 0) = r 5K (t = 0). 

First we calculate the critical density from the equation v = 0. To this end we take into 
account (|33|) . (l39l) and the formulas of Appendix D. As a result, we obtain the dependence 
of the dimensionless critical density r\ c (r) = npa 3 /6) on the ionicity X = y 2 (see Fig. 1). 
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FIG. 1: Critical density as a function of X. 



In order to calculate the chemical potential at the critical point we introduce Av 



Av s — vf^ F i where 



4f = -(N)hsWs(O) 



is the mean-field value of the chemical potential Vs{0). Av c is obtained from the condition 
h = 0; taking into account ( 1351) it yields : 

2 



V 



3& 



[\-zY-2z 



where 



(44) 



S n ( Vc ;0)=G n /(N) 



HS 



(45) 



is the nth particle structure factor at the critical density r] c (X) when ki = 0. In Fig. 2 Az/ C 
is displayed as a function of the ionicity for different values of the parameter z. 

Now we calculate r 2 , u, t and t G at r\ = r\ c using (|32j) . (|34|) . (!38l) - (l39j) . (142]) and formulas 
from Appendices B-D. 

The dependence of r 2 on X at different values of the parameter A is plotted in Fig. 3. 
The explicit formula for Tg R is given in Appendix C. The coefficient u and the shift in the 
mean-field critical temperature, to, as functions of X are plotted in Figs. 4 and 5. As is seen, 
quantities r 2 , to and u are increasing functions of X in the whole region under consideration 



and their dependencies of X are at variance with those obtained in 



361 ] for the lattice model. 
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FIG. 2: Af c as a function of X calculated from (|44p for different values of z (rj = ry c ). The inset 
depicts the behavior of Au c close to the origin. 
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FIG. 3: The dependence of r 2 on the ionicity for different A (rj = rj c ). 

Despite this fact, the behavior of the Ginzburg temperature as a function of X calculated 
in this work is qualitatively similar to that found in 36j] (see Figs. 6-8). Moreover, as in 
36J, the behavior of to{Z) becomes nonmonotonic starting with some value of the attraction 
potential range (A in our case). One can see in Fig. 7 that, for A = 2, to first drop off (at 
very small values of the ionicity) then increases slightly and at X ~ 1.23 again starts to 



15 



0.16 



0.12 



0.08 



0.04 




FIG. 4: The dependence of u on the ionicity Z (77 = r/ c ). 




FIG. 5: The reduced shift of the mean-field critical temperature, to, as a function of I at r] = r/ c . 



decrease. In Fig. 8 the ratio of reduced Ginzburg temperatures, tG(I)/tc(0), is shown at 
different values of A. It is worth noting that the non-monotonic behavior of £g(X) becomes 
more pronounced as A increases. 

In Table 2 we compare our results for the ionicity dependence of the Ginzburg temper- 
ature (at A = 1.5) with the results obtained in [36] for the lattice model as well as with 
experimental data for the crossover temperatures t x (data for I and t x are taken from 



361]). The systems (b)-(d) correspond to the same ionic species Bu4NPic within solvents of 



different dielectric constant. As is seen, in this case our results are in good agreement (qual- 
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I 

FIG. 6: The reduced Ginzburg temperature, to, as a function of X at A = 1.5 (rj = rj c ). The inset 
depicts the behavior of tc{X) close to the origin. 
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FIG. 7: The same as in Fig. 6 but at A = 2. 

itative and quantitative) with the experimental findings. The system (d) is Na in NH 3 and, 
of course, might be described by the potential 4> SR (r) with the different attraction range A. 
For instance, for A = 2 we obtain toiX = 6.97) = 0.8 x 10 -2 (see Fig. 7) that correlates with 
the experimental value t x = 0.6 x 10~ 2 
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FIG. 8: The ratio tc(I) Ag(0) as a function of the ionicity at different values of A (rj = r] c 



TABLE II: Experimentally assessed crossover temperature, t x , taken from 36J]: (a) tetra-ra- 
butylammonium picrate (Bu4NPic) in 1-tridecanol; (b) Bu4NPic in 1-dodecanol; (c) Bu4NPic in 
75% 1-dodecanol plus 25% 1,4-butanediol; (d) Na in NH3; (e) tetra-n-pentylammonium bromide 



in water and the reduced Ginzburg temperature, to, found theoretically in 



361 ] and in this work. 



System 


IonicityT 


tx 


*G ([36]) 


to (this work) 


uncharged fluid 





0(1) 


1 


~ 0.09 


(a) 


17.9 


~ 10~ 3 


~ 0.712 


2.7 x 10" 3 


(b) 


16.8 


~ 0.9 x 1(T 2 


~ 0.717 


0.38 x 10~ 2 


(c) 


8.9 


~3x 10" 2 


~ 0.777 


2.5 x 10" 2 


(d) 


6.97 


~ 0.6 x 10" 2 


~ 0.807 


3.7 x 10" 2 


(e) 


~ 1.4 


Oil) 


1 


~ 0.09 



V. SUMMARY 



In this paper we study the reduced Ginzburg temperature as a function of the interplay 
between the short- and long-range interactions. The ionic fluid is modelled as a charge 
asymmetric continuous system that includes additional short-range attractions. The model 
without Coulomb interactions exhibits a gas-liquid critical point belonging to the Ising class 
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of criticality. We derive an effective GLW Hamiltonian for the model whose coefficients have 
the form of an expansion in powers of the ionicity. Using these coefficients we calculate 
a Ginzburg temperature depending on the ionicity. To this end we introduce a specific 
model which consists of charged hard spheres of the same diameter interacting through the 
additional square-well potentials. To study the effect of the interplay between short- and 
long-range interactions we change, besides the ionicity, the range of the square- well potential. 
As a result, we obtain the similar tendency for the reduced Ginzburg temperature as in 



361 ] when the region of the short-range attraction increases i.e., its nonmonotonic charac- 



ter but with different numerical characteristics. However, our results demonstrate a much 
faster decrease of the Ginzburg temperature when the ionicity increases. We found a good 
qualitative and sufficient quantitative agreement with the experimental findings for B^NPic 
in n-alkanols. This confirms the experimental observations that an interplay between the 
solvophobic and Coulomb interactions alters the temperature region of the crossover regime 
i.e., the increase of the ionicity that can be related to the decrease of dielectric constant 
leads to the decrease of the crossover region. We suggest that the quantitative discrepancy 
of the results for to obtained in 



361 ] and in this work could be due to the fact, besides 
the difference in the symmetry of the two models, that the chemical potential (or density) 
dependence of the Hamiltonian coefficients was taken into account explicitly in our case. 

It should be noted that in the approximation considered in this paper only the critical 
chemical potential depends explicitly on the charge magnitude. In order to obtain the charge 
dependence of the other quantities terms of order higher than y 2 should be taken into account 
into the effective Hamiltonian. Finally, we emphasize that the functional representation (JTj)- 
(|H1) allows to consider more complicated models in particular models including charge and 
size asymmetry. 
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VI. APPENDICES 



A. Recurrence formulas for the cumulants Fourier space. 





• ? k n ) 


— G n (k\, k 2 , . . . , k n ) 






• 5 kn) 


= 




m^(h,k 2 ,.. 


■ 5 kn) 


= (i) 2 Pqlc a G n _ 1 (k 1 ,k 2 , . . ■ , 


|k n _i + k n |) 


m^(k h k 2 ,.. 


■ j kn) 


= (if(3 3 / 2 q 3 a c a G n ^(h,k 2 ,.. 


. . , k„_ 2 + k n _i + k n |) 




■ ■ j kn) 


= iif^h^lc^Gn.^h, 


k 2 , \k n - 2 + k n _i + k„|) 



+ (qt c a - 3 [qlc a ] J Gnsih, k 2 , . . . , |k n _ 3 + . . . + k n |) 



where GrJki, k 2 , . . . , k n ) is the Fourier transform of the n-particle truncated correlation 
function [39( of a one-component hard sphere system and summation over repeated indices 
is meant. 

B. The nth-particle structure factors of a one component hard sphere systems in 
the Percus-Yevick approximation 



*<°> = W (46) 



9 (0) - q-^) 7 (l-^-6r/ 2 ) U7] 



(1 - r7) 10 (l - 30 V + Sir? 2 + U0 V 3 + 60^) 

sm ~ (TTW ' m 

(1 - r/) 13 (l - 8577 + 957?7 2 - 1063r/ 3 - 3590r? 4 - 2940r? 5 - 840r? 6 ) 

Ss(0) " (l + 2r^ (49) 

C. Explicit expression for r| R 

Let us write the Ornstein-Zernike equation in the Fourier space 

*<*> = Tzm< (50) 
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where c(k) is the Fourier transform of the Ornstein-Zernike direct correlation function 14 

n 

We have for c(k) in the Percus-Yevick approximation [41] 

pc(k) = —2477 (afc 3 (sin(fc) - A;cos(A;)) + (3k 2 (2ksm(k) - (k 2 - 2) cos(A;) - 2) 

+ ^r/a((4A; 3 - 24A;) sin(Jfe) - (k 4 - 12k 2 + 24) cos(A;) + 24)^ /k 6 , (51) 

where 

_ (l + 2^) 2 rtH + W 

a (i-^) 4 ' p " (i-vr 

From fl50l) and floTj) we get for g 2 

2 _ <5 22 (0) _ (4r/ 6 - 27rf + 84r? 4 - 146r/ 3 + 144r? 2 - 75r? + 16) 
9 ~ 2^(0) ~ °'° 5r/ (l+W 7 ? ' 

Taking into account ( 1381) we have b 2 SR /ws{Q) = 0.1A 2 . As a result, Tg R is as follows 

2 / (4/7 6 — 27?y 5 + 84?7 4 — 146r? 3 + 144/y 2 — 75?y + 16) 2 



7^ = -0.05^ (1 + 2r/)2( ;_ ?/) 4 ^ + ^(14 /o(/y))). (02) 

where t (2/) is given by 



D. Explicit expressions for the integrals used in equations (j31|) - (l35() 

Using = J (dk) we can present 

^A(fc) = - dxx 2 A(x), (53) 
k n ^° 

?w=^r^^ ))a ' (54) 

^A(A;)A( 2 )(fc) = rdxxA(x)(2f 1 (x)+xf 2 (x)), (55) 

^ K{Jy)Hs Jo 

where the following notations are introduced: 

A(x) = sin (x) (x 3 + k* 2 sin (x)) , (56) 
/i(x) = (x 2 (cos (x) x — 3 sin (x))) (x 6 + 2 k* 2 x 3 sin (x) + k* 4 — k* 4 cos 2 (x)) 1 (57) 
/2OE) = —x (x 5 sin (x) + x 2 k* 2 + k* 2 x 2 cos 2 (x) + 6 cos (x) x 4 — 6 k* 2 x sin (x) cos (x) 
-12 x 3 sin (x) + 6 k* 2 — 6 k* 2 cos 2 (x)) (x 9 + 3 k*V sin (x) + 3 k* 4 x 3 
— 3 k* 4 x 3 cos 2 (x) + k* 6 sin (x) — k* 6 sin (x) cos 2 (x)) 1 (58) 

21 



with x = ka and k* = k d o = ^j2Ay 2 r] being the reduced Debye number. 
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